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Abstract 

In this paper, we study how the distribution of net charges in carbon nanotubes can be influenced 
by substrate and external electric fields, using theoretical calculations based on an extension of 
the atomic charge-dipole model. We find that the charge enhancement becomes less significant 
when the tube gets closer to substrate or when the dielectric constant of substrate increases. It is 
demonstrated that net charges can be shifted to one side of tube by longitudinal electric fields and 
the polarity of charges can be locally changed, while transversal fields give much less infiuence on the 
charge enhancement. These properties could be generalized for other metallic or semiconducting 
nano/microwires and tubes. 
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I. INTRODUCTION 



The distribution of electric charges in carbon nanotubes (CNTs) is of interest for their 
future uses in nanoelectromechanical systems (NEMS)^ such as field emission devices,'^ 
sensors^ actuator^ and charge storages Recently, electric force microscopies (EFM) 
have been used to inject and to detect net charges in 

CNTs,™21 electric char ges are found 
to distribute uniformly along CNTs. However, charge accumulation (so-called charge en- 
hancement) at tube ends has been predicted by theoretical studies using density functional 
theorjff^ and classical electrostatics^^ calculations. These predicted properties have been 
then confirmed by Zdrojek et alJ^ in EFM experiments. It was also shown that electric 
charges can be trapped in CNT loops during periods of timd^. Furthermore, charge-induced 
failure^ and structure changes of CNTj^were reported. In one of our previous works, weak 
charge enhancement at the tube ends and its geometry dependence were demonstrated by 
the combination of theoretical calculations and EFM experiments.'^® 

In this paper, we address the issue of the substrate- and electric-field-effects on the charge 
distribution in CNTs, since CNTs are usually deposited on substrate and driven by electric 
fields in a number of nanodevices.'^^'^ It is known that substrate can exert quite strong 
influence on the charge distribution, as discussed recently in Ref .'^ for the case of ions inside 
CNTs. Theoretical calculations have been performed due to the difficulties for accurately 
quantifying this effect in recent experiments. Our calculation results reveal that the charge 
enhancement becomes less significant when substrate gets closer to CNTs, and that the 
enhancement ratio decreases with increasing dielectric constant of substrate. These effects 
on the charge distribution in radical directions are also discussed. Furthermore, we find 
that the charge distribution in CNTs can be significantly modified in external fields. The 
dependence of field strength is demonstrated for both single-walled and multi-walled CNTs 
(SWCNTs and MWCNTs). We note that the properties demonstrated in this paper could 
also apply to semiconducting CNTs, because semiconducting and metallic nanotubes are 
both expected to accept extra charges, from theoreticaP^ and experiment aP^^ points of 
view. 

The charge distribution has been computed using a Gaussian-regularized atomic charge- 
dipole interacting model.'^^'^ It has been developed from the atomic dipole theory of 
ApplequisiP^ and has recently been parameterized for CNTs.'^Zi In this model, each atom 
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is treated as an interacting polarizable point with a free charge, the static equihbrium state 
of charges are determined by minimizing the total electrostatic energy of system. In this 
work, we have extended this model to take the substrate effect into account by including 
surface-induced terms to vacuum electrostatic interacting tensors using the method of mir- 
ror image.^ Compared with classical Coulomb-law-based models in which only the charge 
is considered, this model provides a more accurate description of electrostatic properties of 
CNTs, because not only the net charges, but also the induced dipole, atomic polarizabilities 
and the image charge are taken into account. 

For the outline, our computational model is presented in sec. II. Results for the effects of 
substrate and fields are discussed in sec. Ill and sec. IV, respectively. We draw a conclusion 
in sec. V. The formulation of the surface-induced electrostatic interacting tensors is given in 
Appendix. 



II. COMPUTATIONAL MODEL 



In our calculation, each atom is associated with an electric charge q and an induced dipole 
p as shown in Fig. [T] 



FIG. 1: Schematic of the principle of the charge-dipole model, in which each atom is modeled as 
a net charge q with an induced dipole p. 

The total electrostatic energy JJ^^^'^ for a CNT of atoms can be written as follows: 
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where x is the electron affinity, V and E stand for the external potential and electric field, 
respectively. T and T are the electrostatic interacting tensors. They can be written as T^l^ = 
(l/47r£o) X (l/rij), = -S/rXc-q ^nd = -V,^, ® V^T^^, where r,,, = \r, - r^]. We 
have regularized T and T by a Gaussian distribution in order to avoid divergence problems 
when atoms are too close to each other, as discussed previously in Refs.'^^'^. Note that the 
value of Gaussian charge distribution width R used in this work for free-end atoms is fitted 
to 0.1273 nm (about 1.3 time that of the carbon atom with three chemical bonds) from 
results in a previous study using DFT calculation.^^ 

The equilibrium state of charges and dipoles should correspond to the minimum value 
of JJ'^^'^'^, and hence the derivatives of JJ'^'''^'^ with respect to q and p should be zero. Taking 
this boundary condition as well as total molecular net charge into account with the 
self-energy terms (when i = j), we can obtain the equilibrium configuration of charge and 
dipole by solving linear vectorial equations and + 1 linear scalar equations as follows: 

' N N 

N N 

i=i i=i (2) 

N 

yi = 1,...,N 

where A is a Lagrange multiplier,'^ which is related to the chemical potential of the 
molecule. In case of a CNT close to a substrate, as in EFM experimentsj ^^ * ^^ ^ the distri- 
butions of charges and dipoles are different from those in free space. We have taken this 
boundary condition into account by adding a surface-induced terms T"^ and T"^ to the 
vacuum electrostatic interacting tensors using the method of mirror images.'^ The detailed 
formulation of T™ and can be found in Appendix. Note that the substrate surface is 
assumed to be infinitely plane in this work. 

The structures of CNTs are relaxed by means of energy optimization^* using the con- 
jugated gradient method^ based on the AIREBO (adaptive interatomic reactive empirical 
bond order) potential^. 
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III. INFLUENCE OF SUBSTRATE 



Previous studies show static charge accumulations at tube end d^^ * ^^ * (as shown in the 
insert of Fig.|2]). A well-defined zone of the charge accumulation is required in order to well 
quantify this enhancement effect. In Fig.[2| we can see that the length of charge enhancement 
zone (L*) increases with the tube length (L), if we define this zone as the part where the 
charge density is higher than the average over the whole tube (a"''^). 
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FIG. 2: L stands for the tube length and L* is the length of the enhancement zone in which the 
average charge density is higher than the average a""^^ over the whole tube. The circles present the 
calculated points. Inset: net charge density along a freestanding tube in space. 

Considering that the CNTs used in experiments are usually longer than those used in our 
calculation, we define the length of charge enhancement zone as 20% of L (10%L at each 
tube end). The ratio of charge enhancement!^ is denoted as follows: 



end I ^middle /0\ 

V9 = a jo (3) 

where o'^'^'^ is the average charge density in the enhancement zone (10%L at each tube 
end), and (J IS that at the middle of the tube. We note that Lp is independent of 0"°""^, 
because the local charge densities are proportional to cr™*^ with respect to a constant electric 
potential on the tube surface. 
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FIG. 3: (a) f versus d, for a charged (5, 5) SWCNT (L 27 nm) with o"''^^ = 0.55 x IQ-^e/atom 
and £2 = CO- Inset: net charge density in a nanotube in a semi-infinite space, (b) ip versus 82, for 
the same tube {d = 0.34nm). 

In case of a CNT in a semi- infinite space (e.g. deposited on substrate), net charges will 
be attracted to the tube bottom by opposite image charges appearing on substrate surface, 
as shown in the inset of Fig. [s] (a). This surface effect mainly depends on the tube-surface 
physisorption distanc#^'^ {d) and the dielectric constant of the substrat^^^'^ (52)- Both 
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of them vary with the type of substrate material. To demonstrate the influence of these 
two parameters on charge enhancement, we plot versus d and 62 in Fig. |3] (a) and (b), 
respectively. We can see in Fig. |3] (a) that the charge enhancement becomes less significant 
when the substrate surface gets closer to the tube (when d decreases). From electrostatic 
point of view, the main mechanism of this effect is that the charge distributed area (band) 
in radical direction has been effectively reduced since a part of net charges is attracted to 
the tube bottom, and hence the charge distribution along the tube axis gets closer to that 
along an infinite-long tube, in which the charge distribution is perfectly uniform [if = 1). 
Similar behavior can be contrasted with the situation when 62 increases, as shown by the 
plot in Fig.js] (b). This implies that (p can get higher if one uses small-dielectric constant 
material instead of metal (£2 = 00) in experiments, and that the CNT exhibits the strongest 
charge enhancement in an infinite space (£2 = 0). 
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FIG. 4: Average charge density versus d for y > (the top half part) and y < (the bottom half 
part), respectively, for a close-ended (9,0) SWCNT (L w 12 nm a"-""^ 2.0 x 10~'^e/atom) on a 
metallic surface. Inset: schematic of transversal charge distribution at a tube end. 



The issue about the charge distribution in radical(transversal) direction has rarely been 
discussed in the literature, although it is one of the main mechanisms of discharging phe- 
nomenon observed experimentally.^'^ When a CNT is horizontally deposited upon a sub- 
strate, charges migration is mainly caused in the direction perpendicular to the substrate 
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surface, as a typical mirror effect (see Fig.|4]). Local charge accumulation at the bottom of 
tube ends (open circles) can directly lead to enhanced electron emission.^^ We can also see 
that the top part of the tube {y > 0) even shows opposite electric sign when d < 0.8 nm 
(solid circles). 




FIG. 5: (Color online) Atomic graph of the charge distribution in a (9,0)@(10,10) DWCNT 
(L ~ 10.3 nm, cj"''^ = 0.55 x lO^^e/atom for both two layers) on a metallic surface. Color of the 
atoms is proportional to the local charge density. Dark arrows stand for the local electric fields 
induced by the net charge around the tube ends, their length and color are proportional to the field 
intensities. The maximum atomic charge densities are: a^"'^ = 8 x lO^^e/atom. The maximum 
strength of the local electric fields presented in this figure E^""^ = 4.3 V/nm. 

The issue about charge distribution in MWCNTs is more complicated due to depolariza- 
tion, field screening and electrostatic interactions between layers. To show further details 
about substrate effects, we depict the atomic charge distribution in a double-walled CNT 
(DWCNT) electrically charged in its both inner and outer carbon layers in Fig.js] The mi- 
gration of atomic charges induced by the metallic surface is shown in this figure. We can see 
the enhanced local electric fields around the tube bottom due to the charge enhancement. 
The top of the tube even shows electrically positive since most of net charges (negative) are 
attracted to the tube bottom by the surface images.!^ 
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IV. INFLUENCE OF EXTERNAL FIELDS 

Recent works showed that external electric fields could induce alignment^, 
deformation^^, field emission!^ and conductivity transitiond^ of CNTs. Here we concern 
mainly on how electric fields infiuence the static distribution of net charges in CNTs. 
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FIG. 6: (Color online) Charge distribution in a charged (9,0) SWCNT (L « 12 nm, a"'"'' = 
1.0 X lO^^e/atom), (a) in space, and (b) in an axial uniform external electric field i?^^* = 0.05 
V/nm. Color of atoms is proportional to charge density. Dark arrows stand for local electric fields 
around tube ends, their length and color are proportional to field intensity. 

The charge distribution of a SWCNT in free space is compared with that in an external 
electric field E^^^ in Fig.[6j As expected, net charges are shifted to one side, around which 
local electric fields (fields induced by net charges and dipoles + external fields) are enhanced. 
The magnitude of this polarization effect is roughly proportional to the external field inten- 
sity E and the tube length this implies that, for CNTs used in experiments (usually 
L ~ /im), E^^^ can be hundreds times weaker for producing a similar effects as those shown 
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in Fig.|6| Moreover, we note that E"^^^ used in this work is about two magnitudes weaker 
than that can lead to field emission from our short CNTs.'^ 
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FIG. 7: Charge profile along a charged (5,5) SWCNT (L « 14 nm, a"""^ = 0.9 x IQ-^e/atom). (a) 
In longitudinal (along the tube axis) electric fields E^^*. Each point is calculated as the average 
of 5%L. (b) In transversal (perpendicular to the tube axis) electric fields i*;*™"*. Each point is 
calculated as the average of 10%L. 

To achieve a quantitative comparison, we plot in Fig. [7] charge distribution along a 
SWCNT in an external field E'^^*. It can be seen that the typical [/-like distributioiJ^ 
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in vacuum (solid circles) can be significantly modified by the axial external electric field. 
On the other hand, the infiuence of transversal electric fields is expected to be weak, due 
to the strong anisotropy of the polarizabilities of CNTs. We can see in Fig.|7| (b) that, even 
with very strong field intensities in the order of V/nm, the charge profile does not change 
a lot. In fact, the transversal field mainly infiuences the charge distribution in non- axial 
direction. Moreover, it needs to mention that the average charge density depends on the 
value of unit length taken in the calculation, e.g. the value of charge density represented 
by the solid circles in Fig.[7] (b) is lower than that in Fig.[7] (a), because it is calculated as 
average on every 10%L, instead of that on every 5%L. 
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FIG. 8: top. Charge profile along the tube axis, for the inner tube of a DWCNT (9,0)@(10,10) 
and for a SWCNT (9,0) with the same size. 

For MWCNTs, it has been reported that electric screening plays an important role for 
the field effects.^Sl To demonstrate its infiuence on the charge enhancement, we compare the 
charge distribution in an inner layer of a DWCNT with that of a SWCNT with the same size 
in Fig.jSj In this comparison, it can be seen that the effect of external fields is much weaker 
on the charge distribution in the inner layer of the DWCNT, due to electrostatic screening. 
Furthermore, by comparing the longitudinal charge distribution with the SWCNT, we have 
found that the charge enhancement is lower in the DWCNT, due to the electric repulsive 
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interaction between the two carbon layers. It is a typical " 1 + 1 < 2 " effect .I^SEZI 



V. CONCLUSION 

In summery, influences of substrate and external electric fields on electric charges in 
CNTs have been investigated by extending the charge-dipole polarization model. The results 
obtained are relevant for a better understanding of the distribution and stability of electric 
charges in CNTs in possible experimental situations (e.g. CNTs deposited on a solid surface). 
Local charge enhancement at tube ends is studied as a particular effect. Our results reveal 
that the charge enhancement becomes less significant when the substrate-tube separation 
decreases or when the dielectric constant of substrate increases. Charge delocalization in 
radical direction has been observed as a typical mirror effect in presence of substrate or 
transversal electric fields. Longitudinal external electric fields have been found to have much 
more influence on the charge enhancement than the transversal ones with same intensities. 
Electric screening in MWCNTs is found to influence charge profile in MWCNTs, especially 
in presence of electric fields. In general, these above conclusions could also qualitatively 
apply to other nanowires and tubes, from electrostatic point of view. 
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VII. APPENDIX: SURFACE-INDUCED TERMS OF ELECTROSTATIC INTER- 
ACTION TENSORS 

To take substrate effects into account, we have extended the charge-dipole model of 
Mayer by adding surface-induced terms {T^, and T^) into the vacuum electrostatic 
interaction tensors (Tq, Ti and T2), respectively, using the method of mirror imaged In 
this method, the electric potential on an arbitrary point (x, y, z) induced by the mirror 
image of a point charge qi embedded in a semi-infinite medium ei close to another medium 
52 (see Fig. [9]) can be written as follows: 
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where r.^ is the coordinate of the mirror image g[" and k = (^2 — ^i)/ [47r£:o(£^2 + ^i)] is an 
electrostatic constant. stands for the 0th order interaction tensor for the mirror image. 
It is the Green's function for the vectorial variable Laplace equation. 
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FIG. 9: Schematic of the charge image. 

For our system shown in this Fig. [9] using Cartesian coordinate, the 0th order mirror- 
charge interaction tensor can be written as: 



n 



-k 



(5) 



where 6x = x — Xi, 6y = y — i/i, ( = z + Zi. 

The interaction tensors of the 1st (charge-dipole) and the 2nd (dipole-dipole) order for 
the image charges can be derived from that of the 0th order (charge-charge): 
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We note that the vacuum interaction tensors used in present study are regularized by a 
normal distribution, in order to avoid divergence problems with point charges when atoms 
get too close to each other.l^ However, it is not necessary to regularize the surface-induced 
terms since the distance between the net charge and its images is generally large enough. 
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